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River networks serve as a paradigmatic example of all branching networks. Essential to under- 
standing the overall structure of river networks is a knowledge of their detailed architecture. Here we 
show that sub-branches are distributed exponentially in size and that they are randomly distribut- 
ed in space, thereby completely characterizing the most basic level of river network description. 
Specifically, an averaged view of network architecture is first provided by a proposed self-similarity 
statement about the scaling of drainage density, a local measure of stream concentration. This scal- 
ing of drainage density is shown to imply Tokunaga's law, a description of the scaling of side branch 
abundance along a given stream, as well as a scaling law for stream lengths. This establishes the 
scaling of the length scale associated with drainage density as the basic signature of self-similarity 
in river networks. We then consider fluctuations in drainage density and consequently the num- 
bers of side branches. Data is analyzed for the Mississippi River basin and a model of random 
directed networks. Numbers of side streams are found to follow exponential distributions as are 
stream lengths and inter-tributary distances along streams. Finally, we derive the joint variation 
of side stream abundance with stream length, affording a full description of fluctuations in network 
structure. Fluctuations in side stream numbers are shown to be a direct result of fluctuations in 
stream lengths. This is the last paper in a series of three on the geometry of river networks. 

PACS numbers: 64.60.Ht, 92.40.Fb, 92.40.Gc, 68.70. +w 



I. INTRODUCTION 

This is the last paper in a series of three on the 
geometry of river networks. In the first [l| we exam- 
ine in detail the description of river networks by scaling 
laws |g, ^, g] and the evidence for universality. Addi- 
tional introductory remarks concerning the motivation of 
the overall work are to found in this first paper. In the 
second article we address distributions of the basic 
components of river networks, stream segments and sub- 
networks. Here, we provide an analysis complementary 
to the work of the second paper by establishing a descrip- 
tion of how river network components fit together. As 
before, we are motivated by the premise that while rela- 
tionships of mean quantities are primary in any investi- 
gation, the behavior of higher order moments potentially 
and often do encode significant information. 

Our purpose then is to investigate the distributions 
of quantities which describe the architecture of river 
networks. The goal is to quantify these distributions 
and, where this is not possible, to quantify fluctuations. 
In particular, we center our attention on Tokunaga's 
law 0, ||, |jj which is a statement about network architec- 
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ture describing the tributary structure of streams. Since 
Tokunaga's law can be seen as the main part of a plat- 
form from which all other river network scaling laws fol- 
low p| , it is an obvious starting point for the investigation 
of fluctuations in river network structure. We use data 
from the Schediegger model of random networks Ji"o| and 
the Mississippi river. We find the distributions obtained 
from these two disparate sources agree very well in form. 
We are able to write down scaling forms of all distribu- 
tions studied. We observe a number of distributions to be 
exponential, therefore requiring only one parameter for 
their description. As a result, we introduce a dimension- 
less scale £t, finding it to be sufficient to describe the fluc- 
tuations present in Tokunaga's law and thus potentially 
all river network scaling laws. Significantly, we observe 
the spatial distribution of stream segments to be random 
implying we have reached the most basic description of 
network architecture. 

Tokunaga's law is also intimately connected with 
drainage density, p, a quantity which will be used 
throughout the paper. Drainage density is a measure 
of stream concentration or, equivalently, how a network 
fills space. We explore this connection in detail, showing 
how simple assumptions regarding drainage density lead 
to Tokunaga's law. 

The paper is structured as follows. We first outline 
Horton-Strahler stream ordering which provides the nec- 
essary descriptive taxonomy for river network architec- 
ture. We then define Tokunaga's law and introduce a 
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scaling law for a specific form of drainage density. We 
briefly describe Horton's laws for stream number and 
length and some simple variations. (Both stream order- 
ing and Horton's laws are covered in more detail in p]). 
We show that the scaling law of drainage density may be 
taken as an assumption from which all other scaling laws 
follow. We also briefly consider the variation of basin 
shapes (basin allometry) in the context of directedness. 
This brings us to the focal point of the paper, the identi- 
fication of a statistical generalization of Tokunaga's law. 
We first examine distributions of numbers of tributaries 
(side streams) and compare these with distributions of 
stream segment lengths. We observe both distributions 
to be exponential leading to the notion that stream seg- 
ments are distributed randomly throughout a network. 
The presence of exponential distributions also leads to 
the introduction of the characteristic number £ t and the 
single length-scale £;(„) cx £t. We then study the variation 
of tributary spacing along streams so as to understand 
fluctuations in drainage density and again find the sig- 
nature of randomness. This leads us to develop a joint 
probability distribution connecting the length of a stream 
with the frequency of its side streams. 



II. DEFINITIONS 
A. Stream ordering 

Horton-Strahler stream ordering jy], |l2) breaks a river 
network down into a set of stream segments. The method 
can be thought of as an iterative pruning. First, we 
define a source stream as the stream section that runs 
from a channel head to the first junction with another 
stream. These source streams are classified as the first- 
order stream segments of the network. Next, we remove 
all source streams and identify the new source streams of 
the remaining network. These are the network's second- 
order stream segments. The process is repeated until one 
stream segment is left of order fi. The order of the basin 
is then defined to be fi (we will use the words basin and 
network interchangeably) . 

In discussing network architecture, we will speak of 
side streams and absorbing streams. A side stream is any 
stream that joins into a stream of higher order, the latter 
being the absorbing stream. We will denote the orders of 
absorbing and side streams by fi and v but when referring 
to an isolated stream or streams where their relative rank 
is ambiguous, we will write stream order as to, subscripted 
as seems appealing. 

Central to our investigation of network architecture 
is stream segment length. As in [||, we denote this 
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length by ' for a stream segment of order u>. We will 
also introduce a number of closely related lengths which 
describe distances between side streams. When referring 
to streams throughout we will specifically mean stream 
segments of a particular order unless otherwise indicated. 
This is to avoid confusion with the natural definition of a 



stream which is the path from a point on a network mov- 
ing upstream to the most distant source. For an order uj 
basin, we denote this main stream length by Z w 

Note also that we consider river networks in planform, 
i.e., as networks projected onto the horizontal (or gravi- 
tationally flat) plane. This simplification poses no great 
concern for the analysis of large scale networks such as 
the Mississippi but must be considered in the context of 
drainage basins with significant relief. 



B. Tokunaga's law 

Defining a stream ordering on a network allows for a 
number of well-defined measures of connectivity, stream 
lengths and drainage areas. Around a decade after the 
Strahler-improved stream ordering of Horton appeared, 
Tokunaga introduced the idea of measuring side stream 
statistics [JtJ ^, ^| . This technique arguably provides the 
most useful measurement based on stream ordering but 
has only recently received much attention ||, |l4|, [l5| . 
The idea is simply, for a given network, to count the 
average number of order v side streams entering an order 
[i absorbing stream. This gives (T M ,„), a set of double- 
indexed parameters for a basin. Note that fi > fi > 
v > 1, so we can view the Tokunaga ratios as a lower 
triangular matrix. An example for the Mississippi river 



is shown in Table 
pictorially in Figure |1 
graph. 



The same data is represented 
in what we refer to as a Tokunaga 
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TABLE I: Tokunaga ratios for the Mississippi River. The 
row indices are the absorbing stream orders while the columns 
correspond to side stream orders. Each entry is the average 
number of order v side streams per order [s, absorbing stream. 



Tokunaga made several key observations about these 
side stream ratios. The first is that because of the self- 
similar nature of river networks, the (T^^) should not 
depend absolutely on either of fi or v but only on the 
relative difference, i.e., k = /i — v . The second is that in 
changing the value of k = fi — v, the (T^) must them- 
selves change by a systematic ratio. These statements 
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streams, i.e., the characteristic scale beyond which ero- 
s^ep cannojt 1 rorc finely dissect the landscape . In 
principle" ara nage density may vary from landscape to 
lfQdscg^ic anc also throughout a single region. Below, we 

j^Ltuwi this observation about Tokunaga's law around 
to show that e 11 river network scaling laws may be derived 
frVm a4 expa ided notion of drainage density. 
^v ei^thoup h the number of side streams entering any 
absorbing str ;am must of course be an integer, Tokuna- 
gE^s ra6os art under no similar obligation since they are 

-a^rag^B. Ne^ ertheless, Tokunaga's law provides a good 
sense of the structure of a network albeit at a level of 
avfcrag&. One of our main objectives here is to go further 

_a^LcoQ)sider fluctuations about and the full distributions 

a t rird important observation of Tokunaga is 



that two of Horton's laws follow from Tokunaga's law, 
which we next discuss. 



Horton's laws 



FIG. 1: A Tokunaga graph for the Mississippi River. The 
values are given in Table |. Each point represents a Tokunaga 
ratio (T Ml „). The solid lines follow variations in the order 
of the absorbing stream p while the dotted lines follow unit 
increments in both p and v, the order of side streams. In 
comparison, the Tokunaga graph of an exactly self-similar 
network would have points evenly spaced at In T\ + (/i + v — 
1) In Rt where 1 < u < p = 2, 3, . . . , fi, i.e., all lines in the 
plot would be straight and uniformly spaced, with the dotted 
lines being horizontal. The nature of deviations in scalings 
laws for river networks is addressed in EJ. 



lead to Tokunaga's law: 



(T M , V ) = (T k ) = (Ti) (R T ) 
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(1) 



Thus, only two parameters are necessary to characterize 
the set of T^y. T\ and Rt- 

The parameter T\ > is the average number of side 
streams of one order lower than the absorbing stream, 
typically on the order of 1.0-1.5. Since these side streams 
of one order less are the dominant side streams of the 
basin, their number estimates the basin's breadth. In 
general, larger values of (Ti) correspond to wider basins 
while smaller values are in keeping with basins with rel- 
atively thinner profiles. 

The ratio Rt > 1 measures how the density of side 
streams of decreasing order increases. It is a measure of 
changing length scales and has a simple interpretation 
with respect to Horton's laws which we describe below. 
Thus, already inherent in Tokunaga's law is a generaliza- 
tion of drainage density p. the usual definition of which is 
given as follows. For a given region of landscape with area 
A with streams totalling in length L, p = L/A and has 
the dimensions of an inverse length scale |Tlj . One may 
think of p as the inverse of the typical distance between 



We review Horton's laws j|, O, Q and then show 
how self-similarity and drainage density lead to Tokuna- 
ga's law, Horton's laws and hence all other river network 
scaling laws. 

The relevant quantities for Horton's relations are n u , 
the number of order ui streams, and {IJ), the average 
main stream length (as opposed to stream segment length 

(lu)) of order oj basins. The laws are simply that the 
ratio of these quantities from order to order remain con- 
stant: 



1 / R n and 
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Ri, 
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for u> > 1. Note the definitions are chosen so that all 
ratios are greater than unity. The number of streams 
decreases with order while all areas and lengths grow. 

A similar law for basin areas (llj IH states that 
(du+i) I (o-u) = Ra where (a u ) is the average drainage 
area of an order u> basin. However, with the assump- 
tion of uniform drainage density it can be shown that 
R n = R a [Q so we are left with the two independent 
Horton laws of equation (^) . 

As in Q, we consider another Horton-like law for 
stream segment lengths: 



= R 



.jC). 



(3) 



As we will show, the form of the distribution of the vari- 
able T Mj „ is a direct consequence of the distribution of 
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Tokunaga showed that Horton's laws of stream num- 
ber and stream length follow from what we have called 
Tokunaga's law, equation ([l]). For example, the solu- 
tion of a difference equation relating the n u and the Tk 
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leads to the result R n = At + [A T — 2Rt] where 
At = (2 + Rt + Ti)/2 for £1 — oo and a more com- 
plicated expression is obtained for finite fl pj, |s|, p~5| . 
In keeping with our previous remarks on Ti , this expres- 
sion for R n shows that an increase in T\ will increase R n 
which, since R a = R n , corresponds to a network where 
basins tend to be relatively broader. Our considerations 
will expand significantly on this connection between the 
network descriptions of Horton and Tokunaga. 

III. THE IMPLICATIONS OF A SCALING LAW 
FOR DRAINAGE DENSITY 

We now introduce a law for drainage density based on 
stream ordering. We write for the number of side 
streams of order v per unit length of order p absorbing 
stream. We expect these densities to be independent of 
the order of the absorbing stream and so we will gener- 
ally use p v . The typical length separating order v side 
streams is then l/p„. Assuming self-similarity of river 
networks, we must have 

Pu+l/pu = l/Rp (4) 

where R p > 1 independent of v. 

All river network scaling laws in the planform may 
be seen to follow from this relationship. Consider an 
absorbing stream of order p. Self-similarity immediately 
demands that the number of side streams of order p, — 1 
must be statistically independent of p. This number is 
of course (Ti). Therefore, the typical length of an order 
p absorbing stream must be 

(lj? ) ) = {T 1 )/p^ 1 . (5) 

Using equation (Q) to replace p^-i in the above equation, 
we find 

Tx/pn-i = RpTi/p^-2- (6) 

Thus, T 2 = R p Tx and, in general T k = (R p ) k - 1 T 1 . This 
is Tokunaga's law and we therefore have 

R P = R T . (7) 
Equation ([I|) and equation (||) also give 

(^ s) ) = T1M.-1 = RpTi/Pv-2 = R P (l^i) • (8) 

On comparison with equation (^) , we see that the above 
is our Hortonian law of stream segment lengths and that 

Rp = Ri(»)- (9) 

As i?;( S ) is the basic length-scale ratio in the problem, 
we rewrite equation (Q), our Hortonian law of drainage 
density, as 

Pv+i/p„ = l/Rns). (10) 

The above statement becomes our definition of the self- 
similarity of drainage density. 



IV. BASIN ALLOMETRY 

Given that we have suggested the need for only a sin- 
gle relevant length ratio, we must remark here on basin 
allomctry. Allomctry refers to the relative growth or scal- 
ing of a shape's dimensions and was originally introduced 
in the context of biology [l7| . A growth or change being 
allometric usually implies it is not self-similar. A long- 
standing issue in the study of river networks has been 
whether or not basins are allometric |, |, [l§]. 

Consider two basins described by (£i,Wi) and 
(T2,W / 2) within the same system where Li is a char- 
acteristic longitudinal basin length and Wi a charac- 
teristic width. The basins being allometric means that 
(Wt/Wz) = (Ti/T 2 ) H where H < 1. Thus, two length 
ratios are needed to describe the allometry of basins. If 
we consider basins defined by stream ordering then we 
have the Horton-like ratios Rt and Rw = Rt- Now, 
when reseating an entire basin, streams roughly aligned 
with a basin's length will rescale with the factor Rt 
and those perpendicular to the basin's axis will rescale 
differently with Rw- This creates a conundrum: how 
can basins be allometric (Rt ^ Rw) an d yet individual 
streams be self-similar (Rt — Rw) as implied by Hor- 
ton's laws? 

We contend the answer is that allometry must be 
restricted to directed networks and that self-similarity of 
basins must hold for non-directed networks. This is in 
agreement with Colaiori et al. |l9| who also distinguish 
between self-similar and allometric river basins although 
we stress here the qualification of directedness. Direct- 
ed networks have a global direction of flow in which the 
direction of each individual stream flow has a positive 
component. A basic example is the random model of 
Scheidegger JTTJ] which we describe below. For a direct- 
ed network, Rt = Ri, and the reseating of basin sizes 
matches up with the rescaling of stream lengths regard- 
less of how the basin's width rescales since all streams 
are on average aligned with the global direction of flow. 
Hence, our premise that streams rescale in a self-similar 
way is general enough to deal with systems whose basins 
rescale in an allometric fashion. 

In considering the allometry of basins, we must also 
address the additional possibility that individual stream 
lengths may scale non-trivially with basin length. In this 
case, the main stream length I would vary with the lon- 
gitudinal basin length L as I cx L d . This is typically a 
weak dependence with 1.0 < d < 1.15 || |2FJ|. Note that 
Horton's laws still apply in this case. The exponent d 
plays a part in determining whether or not a basin scales 
allometrically. The exponent H introduced in the discus- 
sion of basin allometry can be found in terms of Horton's 
ratios (or equivalently Tokunaga's parameters) and d as 
H = d\nR n /\nRi - 1 @. 

Thus, for a directed network d = 1 and H < 1 (e.g., 
Scheidegger |I(J) whereas for undirected, self-similar net- 
works H = 1 and d > 1 (e.g., random undirected net- 
works pl|, E2| ) . River networks are in practice often nei- 
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ther fully directed or undirected. Scaling laws observed 
in such cases will show deviations from pure scaling that 
may well be gradual and difficult to detect H. 

V. TOKUNAGA DISTRIBUTIONS 

The laws of Tokunaga and Horton relate averages of 
quantities. In the remainder of this paper, we investigate 
the underlying distributions from which these averages 
are made. We are able to find general scaling forms of a 
number of distributions and in many cases also identify 
the basic form of the relevant scaling function. 




FIG. 2: Scheidegger's random directed networks. Sites are 
arranged on a triangular lattice and stream flow is directed 
down the page. At each site, the flow direction is randomly 
chosen to be in one of the directions shown on the left. The 
dashed box indicates the area "drained" by the local site. 



To aid and motivate our investigations, we examine, 
as we have done in both ]lj] and ||, a simple model of 
directed random networks that was first introduced by 
Scheidegger fltj|| Since we make much of use this model in 
the present work, we provide a self-contained discussion. 
Consider the triangular lattice of sites oriented as in Fig- 
ure |^. At each site of the lattice a stream flow direction 
is randomly chosen between the two possible diagonal 
directions shown. It is therefore trivial to generate the 
model on a large scale, allowing for a thorough investi- 
gation of its river network statistics. The small, tilted 
box with a dashed boundary represents the area drained 
by the enclosed site. As with many discrete-space mod- 
els, the details of the underlying lattice are unimportant. 
On a square lattice, the model's streams would have three 
choices of flow, two diagonals and straight down the page. 
However, the choice of a triangular lattice does simplify 
implementation and calculation of statistics. For exam- 
ple, only one tributary can exist at each site along a 
stream and stream paths and basin boundaries are pre- 
cisely those of the usual discrete-space random walk [£3| . 

Since random walks are well understood, the exponents 
of many river network scaling laws are exactly known 
for the Scheidegger model |24|, ^5| |2(| p7| and analogies 
may also be drawn with the Abelian sandpile model p8[ . 
For example, a basin's boundaries being random walks 



means that a basin of length L will typically have a width 
W oc L 1 / 2 which gives H = 1/2. Since the network is 
directed, stream length is the same as basin length, I = L, 
so we trivially have d = 1. Basin area a is estimated by 
WL cx L 3 / 2 cx Z 3 / 2 so / cx a 2 / 3 giving Hack's law with an 
exponent of 2/3 

Nevertheless, the Tokunaga parameters and the Hor- 
ton ratios are not known analytically. Estimates from 
previous work f§ find T x ~ 1.35, Ri = Rt — 3.00 and 
R n ~ 5.20. Data for the present analysis was obtained on 
L = 10 4 by W = 3x 10 3 lattices with periodic boundaries. 
Given the self-averaging present in any single instance on 
these networks, ensembles of 10 were deemed sufficient. 

We first examine the distributions of Tokunaga ratios 
Tp v and observe a strong link to the underlying distribu- 
tion of ijf '. Both are well described by exponential dis- 
tributions. To understand this link, we next consider the 
distances between neighboring side streams of like order. 
This provides a measure of fluctuations in drainage den- 
sity and again, exponential distributions appear. We are 
then in a position to develop theory for the joint probabil- 
ity distribution between the Tokunaga ratios and stream 
segment lengths and, as a result, the distribution for the 
quantity T^^/l^ and its inverse. In the limit of large fi, 

(s) 

the T Mj „ are effectively proportional to and all fluc- 
tuations of the former exactly follow those of the latter. 

All investigations are initially carried out for the Schei- 
degger model where we may generate statistics of ever- 
improving quality. We find the same forms for all distri- 
butions for the Mississippi data (and for other river net- 
works not presented here) and provide some pertinent 
examples. Perhaps the most significant benefit of the 
simple Scheidegger model is its ability to provide clean 
distributions whose form we can then search for in real 
data. 

Figure || shows the distribution of the number of order 
v = 2 side streams entering an order \i = 6 absorbing 
stream for the Scheidegger model. At first, it may seem 
surprising that this is not a single-peaked distribution 
centered around {T^. v ) dying off for small and large val- 
ues of T^ v . 

The distribution of in Figure || is clearly well 
described by an exponential distribution. This can also 
be seen upon inspection of Figures [|(a) and ||(b) . Fig- 
ure ^(a) shows normalized distributions of T M ,„ for v = 2 
and varying absorbing stream order n = 4, 5 and 6. 
These distributions (plus the one for absorbing stream 
order n — 7) are rescaled and presented in Figure ||(b). 
The single form thus obtained suggests a scaling form of 
the T Mj „ distribution is given by 

P(T,. V ) = {R m )-^F [T^tow)-"] . (11) 

where F is an exponential scaling function. However, 
this only accounts for variations in /z, the order of the 
absorbing stream. 

Figures |5|(a) and ||(b) show that a similar rescaling of 
the distributions may be effected when v is varied. In this 
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FIG. 3: An example of a generalized Tokunaga distribu- 
tion for the Scheidegger model. The Tokunaga ratio T^^ is 
the number of side streams of order u entering an absorbing 
stream of order fj,. For this particular example n = 6 and 
v — 2. The form is exponential and is a result of variations 
in stream segment length rather than significant fluctuations 
in side stream density. 



case, the data is for the Mississippi. The rescaling is now 
by i?;( a ) rather than R7(^ and equation (|ll]) is improved 
to give 

P(T M ,„) = (R lW r-"- 1 P r [T^/iR^y-"- 1 ] . (12) 

The function Pt is a normalized exponential distribution 
independent of /x and v, 



Pt{z) = -e^* 
Si 



(13) 



where £t is the characteristic number of side streams of 
one order lower than the absorbing stream, i.e., £t = 
(Ti). For the Mississippi, we observe £t ~ 1.1 whereas 
for the Scheidegger model, £t ~ 1.35. As expected, the 
Tokunaga distribution is dependent only on k = fi — v so 
we can write 

P(T k ) = (R^-'Pt [T^/(R l(s) ) k -i] . (14) 

with Pt as above. 



VI. DISTRIBUTIONS OF STREAM SEGMENT 
LENGTHS AND RANDOMNESS 

As we have suggested, the distributions of the Tokuna- 
ga ratios depend strongly on the distributions of stream 
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FIG. 4: Distributions for Tokunaga ratios for varying orders 
of absorbing stream and fixed side stream order of v = 2 
for the Scheidegger network. In (a), examples of T Ml „ dis- 
tributions for absorbing stream order = 4 (circles), fi = 5 
(squares) and fj, — 6 (triangles). In (b), these distributions, as 
well as the fi = 7 case, are rescaled according to equation (|l 
The resulting "data collapse" gives a single distribution. For 
the Scheidegger model, i?,( s ) — 3.00. 



segment lengths. Figure ^| is the indicates why this is so. 
The form of the underlying distribution is itself exponen- 
tial. We have already examined this fact extensively in || 
and here we develop its relationship with the Tokunaga 
distributions. 

Figures |(a) and |(b) show that the distributions of 
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Ifi can be rescaled in the same way as the Tokunaga 
distributions. Thus, we write the distribution for stream 



segment lengths as M 

P(l^) = (R m )-" +1 P m fc s) /0R/( 3 > 



-ii+i 



(15) 



As for Pt, the function P;( B ) 
distribution 



is a normalized exponential 



P lls) ( z ) = -L e - Z ^K, 

SiM 



(16) 



where, in a strictly self-similar network, £ ; ( a ) is the charac- 
teristic length of first-order stream segments, i.e., £;( S ) = 
Hi^}. (Note that in |j) we use £ for £ ;(s) for ease of 
notation). We qualify this by requiring the network to 
be exactly self-similar because in most models and all 
real networks this is certainly not the case. As should be 
expected, there are deviations from scaling for the largest 
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FIG. 5: Tokunaga distributions for varying side stream 
orders for the Mississippi river basin. In both (a) and (b) , the 
absorbing stream order is fi = 5 and the side stream orders 
are v = 2 (circles), v = 3 (squares) and v — 4 (triangles). 
The raw distributions are shown in (a). In (b) the distribu- 
tions are rescaled as per equation (|12|). For the Mississippi, 
the ratio is estimated to be R t (a) — 2.40 M. 
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FIG. 6: Stream segment length distributions for varying 
stream order for the Scheidegger model. Lengths are in units 
the lattice spacing. Shown in (a) are raw distributions for lo — 
4 (circles), lu — 5 (squares) and lo = 6. The linear forms on the 
semilogarithmic axes indication these distributions are well 
approximated by exponentials |J. In (b), the distributions in 
(a) plus the distribution for Ij 1 (diamonds) are rescaled using 
equation (119). 



Oh 

3 



and smallest orders. Therefore, £j( 3 > is the characteristic 
size of a first-order stream as determined by scaling down 
the average lengths of those higher order streams that are 
in the self-similar structure of the network. It is thus in 
general different from ). 

We therefore see that the distributions of and ijf* 
are both exponential in form. Variations in largely 
govern the possible values of the T^ v . However, T^ tV 

(s) 

is still only proportional to on average and later on 
we will explore the joint distribution from which these 
individual exponentials arise. 

The connection between the characteristic number £ t 
and the length-scale follows from equations (||), @, 
and @: 

Zt = PiRi(*)t;m- (17) 

This presumes exact scaling of drainage densities and in 
the case where this is not so, p\ would be chosen so that 
(-Rj(s) ) v ~ x p\ most closely approximates the higher order 

Pv 

We come to an important interpretation of the expo- 
nential distribution as a composition of independent 
probabilities. Consider the example of stream segment 
lengths. We write as the probability that a stream 



segment of order /i meets with (and thereby terminates 
at) a stream of order at least fj,. For simplicity, we assume 
only one side stream or none may join a stream at any 
site. We also take the lattice spacing a to be unity so that 
stream lengths are integers and therefore equate with the 
number of links between sites along a stream. For a =/= 1, 
derivations similar to below will apply with replaced 

by [ijy /a], where [•] denotes rounding to the nearest 
integer. Note that extra complications arise when the 
distances between neighboring sites are not uniform. 

Consider a single instance of an order /x stream seg- 
ment. The probability of this segment having a length 
is given by 

PG«)=iUi-iv)^ (is) 

where p M is the probability that an order an order fj, 
stream segment terminates on meeting a stream of equal 
or higher order. We can re-express the above equation as 

P(^))^p AI exp{-?( s )ln(l-p,)- 1 }, (19) 

and upon inspection of equations ( |l5| ) and ( |l6| ) we make 
the identification 

(R^r-^m = [-Mi-P»T\ (20) 



8 




50 100 150 200 

FIG. 7: A comparison of inter-tributary length distributions 
for the Scheidegger model. The example here is for order 
H — 6 absorbing streams and order v = 3 side streams. Note 
that no rescaling of the distributions has been performed. The 
three length variables here correspond to x=b, x=i and x=e, 
i.e., /{f,'„ b ' (circles) ify (squares), and (triangles). These 
are the beginning, internal and end distances between enter- 
ing side branches, defined fully in the text. No quantitative 
difference between these three lengths is observed. 



which has the inversion 

p /1 = l- e -V(*,w)'*- 1 € l w. (21) 

For (i sufficiently large such that <C 1, we have the 
simplification 

p^i/^r 1 ^. (22) 

We see that the probabilities satisfy the Horton-like scal- 
ing law 

pVpm-i = V-Rjm- ( 23 ) 

Thus, we begin to see the element of randomness in our 
expanded description of network architecture. The ter- 
mination of a stream segment by meeting a larger branch 
is effectively a spatially random process. 

VII. GENERALIZED DRAINAGE DENSITY 

Having observed the similarity of the distributions of 
Tfj, M and , we proceed to examine the exact nature of 
the relationship beteween the two. To do so, we introduce 
three new measures of stream length. These are 1^,'^, 




FIG. 8: An examination of the asymptotic behavior of dis- 
tributions of internal inter-tributary lengths. The data here 
is for the Scheidegger model for the case of fixed side stream 
order v = 2. The plot in (a) shows distributions for absorb- 
ing stream order fi = 3 (circles), fi — 4 (squares), fi = 6 
(triangles), and fi — 8 (diamonds). As fi increases, the distri- 
butions, which are all individually exponential, tend towards 
a fixed exponential distribution. Since lower order stream seg- 
ments have typically smaller lengths, they statistically block 
larger values of £^,'„ , reducing the extent of the distribution 
tails for low /j,. This is further evidenced in (b) which pro- 
vides a plot of (Z^% ), the mean inter-tributary stream length, 
as a function of [i with v — 2. These mean values approach 
{lu'uLz) ~ V/°2 — 10.1 where p2 is the density of second-order 
side streams. 



the distance from the beginning of an order /i absorbing 
stream to the first order v side stream; , the distance 
between any two adjacent internal order v side streams 
along an order [i absorbing stream; and I 8 ' the distance 
from the last order v side stream to the end of an order 
fi absorbing stream. By analysis of these inter-tributary 
lengths, we will be able to discern the distribution of 
side stream location along absorbing streams. This leads 
directly to a more general picture of drainage density 
which we fully expand upon in the following section. 

Figure compares normalized distributions of ifcv , 
l^'v and I 8 ' % for the Scheidegger model. The data is for 
the distance between order v = 3 side streams entering 
order fi = 6 absorbing streams. Once again, the distribu- 
tions are well approximated by exponential distributions. 
Moreover, they are indistinguishable. This indicates, at 
least for the Scheidegger model, that drainage density is 
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independent of relative position of tributaries along an 
absorbing stream. 

We now consider the effect on the distribution of inter- 

fs i) 

nal inter-tributary distances l^v following from varia- 
tions in /z, the order of the absorbing stream. Figure H(a) 

(si). 

provides a comparison of l^u distributions for v = 2 
and fi — 3 through fi = 8. As /i increases, the distri- 
butions tend towards a limiting function. With increas- 
ing /i we are, on average, sampling absorbing streams of 

fs i) 

greater length and the full range of 1^'^' becomes accord- 
ingly more accessible. This approach to a fixed distri- 
bution is reflected in the means of the distributions in 
Figure ||(a). Shown in Figure ||(b), the means (Z^f/iP) for 
v = 2 approach a value of around 10.1. The correspond- 
ing density of second-order streams for the Scheidegger 
model is thus pi = l/il^'vL^) — 0.01. Higher drainage 
densities follow from equation (|l0|). However, since devi- 
ations occur for small v, there will also be an approach 
to uniform scaling to consider with drainage density. 



VIII. JOINT VARIATION OF TOKUNAGA 
RATIOS AND STREAM SEGMENT LENGTH 

We have so far observed that the individual distribu- 

fs) 

tions of the ' and T Mjiy are exponential and that they 
are related via the side-stream density p n u. However, 
this is not an exact relationship. For example, given a 

(s) 

collection of stream segments with a fixed length ' we 
expect to find fluctuations in the corresponding Tokuna- 
ga ratios T^„. 

To investigate this further we now consider the joint 

(s) 

variation of T^ v with from a number of perspectives. 
After discussing the full joint probability distribution 
P(Tfj.^, Ip) we then focus on the quotient v = T^^/ljf^ 
and its reciprocal w = ijf /T^ tV . The latter two quanti- 
ties are measures of drainage density and inter-tributary 
length for an individual absorbing stream. 



A. The joint probability distribution 

We build the joint distribution of P(T Miiy , iff') from 
our conception that stream segments are randomly dis- 
tributed throughout a basin. In equation (Q), we have 

fs) 

the probability of a stream segment terminating after 
steps. We need to incorporate into this form the prob- 
ability that the stream segment also has order v 
side streams. Since we assume placement of these side 
streams to be random, we modify equation (|l^) to find 

P(l<?\T», v ) =P^l~ l )p T ^{l-Pu-p^- T ^-\ 

(24) 

where (?) = n\/k\(n — k)\ is the binomial coefficient and 
p v is the probability of absorbing an order v side stream. 



The extra p v appears in the last factor (1 — p v — p^) 
because this term is the probability that at a particular 
site the stream segment neither terminates nor absorbs 
an order v side stream. Also, it is simple to verify that the 
sum over iffi and T^ v of the probability in equation ( |l9| ) 
returns unity. 

While equation ( |l9| ) does precisely describe the joint 
distribution P(l^ s \ T^, u ), it is somewhat cumbersome to 
work with. We therefore find an analogous form defined 
for continuous rather than discrete variables. We simplify 
our notation by writing p = p v , q = (1 — p v — p^) and 

p = p^. We also replace (l^jT^) by (x,y) where now 
i,j/£l, Note that < y < x — 1 since the number of 
side streams cannot be greater than the number of sites 
within a stream segment. 
Equation ( |l9| ) becomes 



P(x,y)=Np 



T(x) 



r(y + l)T(x-y) 



(25) 



where we have used T(z + 1) = z\ to generalize the 
binomial coefficient. We have included the normaliza- 
tion N to account for the fact that we have moved to 
continuous variables and the resulting probability may 
not be cleanly normalized. Also we must allow that 
N = N(p,p) and we will be able to identify this form 
more fully later on. Using Stirling's approximation p9[ , 
that T(z + 1) ~ V27rz z+1/2 e~ z , we then have 



P(x,y) = Nppyqz-y- 1 - 



1 



(X - lf-3/ 2 



2tt y y+1 / 2 (x -y - l)z-2/-i/2 



~ N'x" 1 / 2 [F(y/x)] x 

where we have absorbed N and all terms involving only 
p and p into the prefactor N' = N'(p,p) = Np/(^/2irq). 
We have also assumed x is large such that x — 1 ~ x and 
l>l/a;~0. 

The function F(v) = F(v;p, q) identified above has the 
form 



F(v) = 



1 



-(l-i,) 



(27) 



where < v < 1 (here and later, the variable v will refer 
to y/x). Note that for fixed x, the conditional probability 
P(y | x) is proportional to [F{y/x)] x . Figure ||(a) shows 
[^(v)]^ for a range of powers x. The basic function has a 
single peak situated near v = p. For increasing x which 

fs) 

corresponds to increasing , the peak becomes sharp- 
er approaching (when normalized) a delta function, i.e., 
]im xroo [F(v)]* = S(v-p). 

Figure |^(b) provides a comparison between data for the 

Scheidegger model and the analytic form of P(lff , T^ J). 
For this example, fj, = 6 and v = 2 which corresponds to 
p ~ 0.10, q ~ 0.90 and p ~ 0.001 (using the results 
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(s) 



v = T // 



FIG. 9: Form of the joint distribution of Tokunaga ratios and 
stream segment lengths. The distribution is given in equa- 
tion ( jir^ ) and is built around the function F(v — T^^/l^) 
given m equation @. Shown in (a) is [F(v)} 1 ^ for = 1, 
10, 100 and 1000. Increasing 1^ corresponds to the focusing 
of the shape. In (b), the distribution P(T M ,„ | ~ 340) is 
compared between theory (smooth curve) and data from the 
Scheidegger model (circles). The Scheidegger model data is 
compiled for a range of values of iji' rescaled as per equa- 
tion (p_h. 



of the previous section). The smooth curve shown is 
the conditional probability P(y | X) for the example val- 



ue of X 



340 following from equation (26). 



From simulations, we obtain a discretized approximation 
to P(lji ,T^ V ). For each fixed x = in the range 

165 < If' < 345, we rescale the data using the following 
derived from equation 



P{X, y) = N'X- 1 / 2 (jV'~ V/ 2 P(x, yj\ 



X/x 



(28) 



All rescaled data is then combined, binned and plotted as 
circles in Figure ^|(b), showing excellent agreement with 
the theoretical curve. 



B. Distributions of side branches per unit stream 
length 

Having obtained the general form of P(lff\ T^ >v ), we 
now delve further into its properties by investigating the 



distributions of the r&tiQbt) = T^/lp and its reciprocal 
w. 

Th&M 



^ 0.5 

'Agcordir 
-cation 



> First, 



We trar 



tern dcstiribcdQlqg^u, v) with^hp relations Q 15 
u 2 = x 2 + i$ = Tand / / v = y/x. 

The inverse relations are x = u/(l + v 2 ) and y — uv/(l + 
v 2 ) and we also have dxdy = xdudv. Equation ( |32| ) leads 
to 



[lantity T^ v /lji' isAhc Qumber of side streams 
per leng' h of a given absoToSig steeam and when aver- 
aged ove : an ensemble oFabsorbingVtreams gives 



Pv 



gly, the reciprocal l^ 1 jT^ v is tlite average 
.side streams of order v. 
_ . we deriWP^/iW) from P(l^T^ v ). 
Pttien consider some intuitive rescalings which will 
us te-0l&:uce the form of the normalization N(p, q) 
We re 1 mte equation (Eq) as 



V) 



-1/2 



exp{— xln [—F(y/x)]} AO 



form (jx, y) to the modified polar coordinate 



(29) 

epa- 
We 

How 



(30) 
sys- 

(31) 



\l + v 



1/2 

.2 I 



■ln[-F(v)} 



(32) 



To find P(v) we integrate out over the radial dimension 



P(v) 



N' 



duP(u, v) 



du 



1 



v 



1/2 
o ; ex P 



= N'(l + v 2 ){\n[-F(v)})- 3/2 
1 + v 2 



dzz 1 ' 2 ^ 



■In [-F{v)\ 



2 = 



= N 



(\n[-F(v)}) 3 / 2 



(33) 



Here, N" = N'T (3/ 2) = and we have used the 

substitution z = u/(l + v 2 ) ln[— F{v)\. 

The distribution for w = ijf /T^.v = 1/v follows sim- 
ply from equation (|33|) and we find 



P(w) = N' 



1 



w 4 (ln[-F(l/w)]) 3 / 2 



(34) 



Figures |l(i|(a) and |Io](b) compare the predicted forms of 
P(v) and P(w) with data from the Scheidegger model. In 
both cases, the data is for order v = 2 side streams being 
absorbed by streams of order \i = 6. Note that both dis- 
tributions show an initially exponential-like decay away 
from a central peak. Moreover, the agreement is excel- 
lent, offering further support to the notion that the spa- 
tial distribution of stream segments is random. 
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FIG. 10: Comparison of theory with measurements of aver- 
age inter-tributary distances for the Scheidegger model. The 
data in both (a) and (b) is for the case of order v = 2 side 
streams and order p = 6 absorbing streams. In (a), the distri- 
bution of v = Tfn^/lfi obtained from the Scheidegger model 
(circles) is compared with the smooth curved predicted in 
equation (^). The same comparison is made for the recipro- 
cal variable w — /T^v, the predicted curve being given in 
equation (0). 



Finally, we quantify how changes in the orders p and v 
affect the width of the distributions by considering some 
natural rescalings. Figure pT|(a) shows binned, normal- 
ized distributions of T^u/ljf' for the Scheidegger model. 
Here, the side stream order is v — 2 and the absorbing 
stream orders range over p = 5 to p = 8. All distribu- 
tions are centered around p 2 — 0.10. 

(s) 

Because the average length of increases by a factor 
i?;( 3 ) with p, the typical number of side streams increases 

(s) 

by the same factor. Since we can decompose as 

l[f ) =l^+l^+--- + l^ + ^, (35) 

where there are T„ „ — 1 instances of l^v\ becomes 

better and better approximated by + 1)(Z^% ). 

Hence, the distribution of T^ v /ljy peaks up around 
P2 as p incr eases, the typical width reducing by a factor 
of Ri(e) for every step in p. Using this observation, 
Figure |ll](b) shows a rescaling of the same distributions 
shown in Figure [Tl](a). The form of this rescaling is 

p(w^ s) ) = (R^r /2 Gi (pwjw - p2]{r^y /2 ) 

(36) 



FIG. 11: Distributions of the quantity T Mi „/i,i for the Schei- 
degger model with p, the order of the absorbing stream, vary- 
ing and the side stream order fixed at v = 2. Given in (a) are 
unrescaled distributions for p — 5 (circles), p = 6 (squares), 
jit = 7 (triangles), and fi = 8 (diamonds). Note that as the 
order of the absorbing stream increases so does its typical 
length. This leads to better averaging and the standard devi- 
ation of the distribution decays as R l 1 . The distributions 
are all centered near the typical density of order v = 2 side 
streams, p2 — 0.10. The rescaled versions of these distribu- 
tions are given in (b) with the details as per equation (prj). 



where the function is similar to the form of P(v) given 
in equation (^3|). The mean drainage density of p2 has 
been subtracted to center the distribution. 

We are able to generalize this scaling form of the distri- 
bution further by taking into account side stream order. 
Figures |l2|(a) and |l2|(b) respectively show the unrescaled 

and rescaled distributions of T^v/ljf with v allowed to 
vary. This particular example taken from the Scheideg- 
ger model is for p — 6 and the range v = 1 to v = 5. 
Since v is now changing, the centers are situated at the 
separate values of the p v . Also, the typical number of side 
streams changes with order v so the widths of the distri- 
butions dilate as for the varying p case by a factor yiiiW . 
Notice that the rescaling works well for v = 2, . . . ,5 but 
not v = 1. As we have noted, deviations from scaling 
from small orders are to be expected. In this case, we 
are led to write down 

P(T^/l^)(R l(B) y»/ 2 G 2 ([T M , v /lU - PuMRii*))-^) 

(37) 

where, again, G^z) is similar in form to P(v). 
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FIG. 12: Distributions of number of side streams per unit 
length for the Scheidegger model with v, the order of side 
streams, varying. For both (a) and (b), the absorbing stream 
order is fj, = 6. Shown in (a) are the unrescaled distribu- 
tions for v — 2 (circles), v = 3 (squares), and v = 4 (tri- 
angles). Note that as v increases, the mean number of side 
streams decreases as do the fluctuations. The distributions in 
(a) together with the distribution for v = 5 (diamonds) are 
shown rescaled in (b) as per equation (p8[). 
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FIG. 13: Tokunaga statistics for the Mississippi river basin. 
The distributions are as per Figure p"2|(a), distributions of 
number of side streams per unit length with v, the order of 
side streams, varying. The absorbing stream order is n — 
7 and the the individual distributions correspond to v = 2 
(circles), v = 3 (squares) and v = 4 (triangles). All lengths 
are measured in meters. Rescalings of the distributions shown 
in (a) along with that for v = 5 (diamonds) are found in (b). 
Reasonable agreement with equation (B8|) is observed. 



We find the same rescalings apply for the Mississip- 
pi data. For example, Figure [l3 (a) shows unrescaled 

distributions of Tp tU /lu for varying v. Figure |l3|(b) 
then shows reasonable agreement with the form of equa- 
tion (|37|). In this case, the Scheidegger model clearly 
affords valuable guidance in our investigations of real 
river networks. The ratio R^s) — 2.40 was calculated 
from an analysis of f<i and l u . The density p2 — 0.0004 
was estimated directly from the distributions oiT^^/l^ 
and means that approximately four second-order streams 
appear every ten kilometers. 

Combining equations ( |38| ) and ( |38| ) , we obtain the com- 
plete scaling form 

P(T M ,„/i W) = {R lW )^ v -WG ([T^/li s) - pA(Ri^) (m 

(38) 

As per G\ and G2, the function G is similar in form to 
P(v). 

The above scaling form makes intuitive sense but is 
not obviously obtained from an inspection of (|33|). We 
therefore examine P(v) by determining the position and 
magnitude of its maximum. Rather than solve P'(v) = 
directly, we find an approximate solution by considering 



the argument of the denominator, — \nF(v), with F(v) 
given in equation j2^). Since the numerator of P(v) is 
1 + v 2 and the maximum occurs for small v this is a 
justifiable step. Setting dF/dv = 0, we thus have 



In- 



1 



'1 



+ lavp = 0, 



(39) 



which gives v m = p/(q + p) = p/(l — p). Note that for 
p < 1, we have v m ~ p. 

Substituting v = v m — p/(l — p into equation (|33|), we 
find 



P(v m ) ~ N"p- 3 / 2 = Np- 112 ^ 12 



(40) 



""Tffeiuming p 2 « 1 and q ~ 1. Returning to the scaling 
forrn of equation (|38|), we see that the p -1 / 2 factor in 
equation ( pp[ ) accounts for the factors of (Rim)^^ 2 since 
p = p~n scales from level to level by the ratio -R;(s). We 
therefore find the other factor (R^b))"^ 2 of equation j3^ ) 
gives N — cp 1 ! 2 where c is a constant. Since p = p u , 
it is the only factor that can provide this variation. We 
thus have found the variation with stream order of the 
normalization N and have fully characterized, P(x,y), 
the continuum approximation of P(lfi , T^J). 
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IX. CONCLUDING REMARKS 

We have extensively investigated river network archi- 
tecture as viewed in planform. We identify the self- 
similarity of a form of drainage density as the essence 
of the average connectivity and structure of networks. 
From previous work in we then understand this to 
be a base from which all river network scaling laws may 
obtained. 

We have extended the description of tributary struc- 
ture provided by Tokunaga's law to find that side stream 
numbers are distributed exponentially. This in turn is 
seen to follow from the fact that the length of stream 
segments are themselves exponentially distributed. We 
interpret this to be consequence of randomness in the 
spatial distribution of stream segments. Furthermore, 
the presence of exponential distributions indicate fluc- 
tuations in variables are significant being on the order 
of mean values. For the example of stream segment 
lengths, we thus identify ^< H ), a single parameter need- 
ed to describe all moments. This is simply related to 
£t, which describes the distributions of Tokunaga ratios. 
The exponential distribution becomes the null hypothe- 
sis for the distributions of these variables to be used in 
the examination of real river networks. 

We are able to discern the finer details of the connec- 
tion between stream segment length and tributary num- 
bers. Analysis of the placement of side streams along a 
stream segment again reveals exponential distributions. 
We are then able to postulate a joint probability dis- 
tribution for stream segment lengths and the Tokunaga 
ratios. The functional form obtained agrees well with 
both model and real network data. By further consider- 
ing distributions of the number of side streams per unit 
length of individual stream segments, we are able to cap- 
ture how variations in the separation of side streams are 



averaged out along higher-order absorbing streams. 

By expanding our knowledge of the underlying dis- 
tributions through empiricism, modeling and theory, we 
obtain a more detailed picture of network structure with 
which to compare real and theoretical networks. We have 
also further shown that the simple random network mod- 
el of Scheidegger has an impressive ability to produce 
statistics whose form may then be observed in nature. 
Indeed, the only distinction between the two is the exact 
value of the scaling exponents and ratios involved since 
all distributions match up in functional form. 

We end with a brief comment on the work of Cui et 
al. |U| who have recently also proposed a stochastic gen- 
eralization of Tokunaga's law. They postulate that the 
underlying distribution for the T^ v is a negative binomial 
distribution. One parameter additional to T\ and Rt, a, 
was introduced to reflect "regional variability," i.e., sta- 
tistical fluctuations in network structure. This is in the 
same spirit as our identification of a single parameter £ t . 
However, our work disagrees on the nature of the under- 
lying distribution of T^v We have consistently observed 
exponential distributions for T^^ in both model and real 
networks. 

In closing, by finding randomness in the spatial distri- 
bution of stream segments, we have arrived at the most 
basic description of river network architecture. Under- 
standing the origin of the exact values of quantities such 
as drainage density remains an open problem. 
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